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Abstract: We define a class of "optimal" coordinate systems by requiring that the 
deviation from an exact Robertson- Walker metric is "as small as possible" within a 
given four dimensional volume. The optimization is performed by minimizing several 
volume integrals which would vanish for an exact Robertson- Walker metric. Covariance 
is automatic. Foliation of space-time is part of the optimization procedure. Only the 
metric is involved in the procedure, no assumptions about the origin of the energy- 
momentum tensor are needed. A scale factor does not show up during the optimization 
process, the optimal scale factor is determined at the end. The general formulation is non 
perturbative. An explicit perturbative treatment is possible. The shifts which lead to the 
optimal coordinates obey Euler-Lagrange equations which are formulated and solved in 
first order of the perturbation. The extension to second order is sketched, but turns out to 
be unnecessary. The only freedom in the choice of coordinates which finally remains are 
the rigid transformations which keep the form of the Robert son- Walker metric intact, i.e. 
translations in space and time, spatial rotations, and spatial scaling. Spatial averaging 
becomes trivial. In first order of the perturbation there is no backreaction. A simplified 
second order treatment results in a very small effect, excluding the possibility to mimic 
dark energy from backreaction. This confirms (as well as contradicts) statements in the 
literature. 
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1 Introduction 



The averaging problem, i.e. the problem of averaging a realistic inhomogeneous metric 
into a smooth one, as well as the fitting problem, the fitting of an "optimal" Robertson- 
Walker (R-W) metric to a realistic inhomogeneous metric, are both non trivial due to 
the freedom of choosing arbitrary coordinates in general relativity. Most papers focus 
on "gauge transformations", where the R-W background metric is given and fixed, 
and only the perturbation /i^j^ is transformed. To fix the background one usually resorts 
to a flow of matter. This is unsatisfactory for two reasons. Firstly one has to make 
rather stringent assumptions concerning the flow, like only one single component and 
absence of rotation. The second point, although quite obvious but nevertheless hard to 
find mentioned in the literature, appears even more drastic: The flow which is used is, of 
course, not the real flow of matter but already some average over an irregular flow. The 
result of this averaging clearly depends on the choice of coordinates. One is faced with 
the bizarre situation that one starts some sophisticated "gauge invariant" averaging 
procedure on the basis of a background obtained from an ambiguous and unspecifled 
averaging. This makes the whole procedure quite dubious. 

There is an extensive literature on both topics which can only be briefly addressed 
here. We refer e.g. to the monograph of Krasihski [T] and the comprehensive review of 
Buchert [2j. A careful analysis of the fltting problem was given by Ellis and Stoeger [5]. 
We only briefly mention some aspects here. The averaging problem was flrst raised by 
Shirokov and Fisher |3] in 1963. The authors suggested to integrate the metric tensor 
over a four dimensional volume with the familiar factor ^—g in the measure. Such an 
expression is, however, not covariant for a tensor due to the freedom of performing local 
transformations. A covariant averaging prescription can be constructed by introducing 
a bivector g^{x, x') of geodesic parallel displacement, as discussed in the appendix of 
This transforms as a vector with respect to coordinate transformations at either x or 
x' and maps a vector Api^x') to Aa{x) = g^{x,x')Ai3{x'), analogously for higher order 
tensors. An averaging with the help of bivectors is also used in the work of Zalaletdinov 
[6J where the emphasis was on the commutativity of averaging and covariant differenti- 
ation. As remarked by Stoeger, Helmi, and Torres |17j the method of using a covariantly 
conserved bivector is not applicable to the metric, because the covariant derivative of 
the metric vanishes. The metric is therefore invariant under this averaging procedure. 
Another popular method due to Bardeen |8] is to work with gauge invariant (in flrst 
order of the perturbation and for static transformations only) quantities. The most gen- 
eral covariant and translation invariant flrst order averaging scheme has been given in 
[9]. But any such an averaging has the principal problem that a plane wave, instead of 
being averaged to zero, will always stay a plane wave, albeit with reduced amplitude: 
exp(ikx) — > / / (x — y) exp(iky)(i^?/ = { / /(z) exp{—ikz)d^z} ■ exp(ikx). 

Instead of attempting an averaging, it appears therefore more promising to deter- 
mine directly an "optimal" approximating smooth metric. Our approach is conceptually 
simple. We flx the coordinate system as far as principally possible, so that no unphysical 
gauge freedom remains. The coordinate system is chosen in such a way that, in a given 
four dimensional volume, the metric is as close to an exact R-W metric as possible. 
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Before going into details one should recall that an exact R-W metric (with A; = 0) keeps 
its form under an eight parameter group of global symmetry transformations: Rigid 
translations in space and time, rigid spatial rotations, and rigid scaling of the spatial co- 
ordinates. (For k = ±1 there is only a seven parameter group, scaling is not allowed.) In 
the case of translations in time, and of scaling in space, the scale function a{t) changes. 

This freedom in the choice of coordinates cannot and need not be fixed. It is inevitably 
connected with the symmetry of the R-W metric. Any covariant fitting procedure will 
necessarily share this freedom of transformations. Therefore a maximal fixing of the 
coordinate system means that the coordinates are fixed up to the above rigid transfor- 
mations, while no further transformations are allowed anymore. The coordinate system 
obtained at the end should fulfill two criteria: 

Covariance: Let two observers A, B describe the same reahstic inhomogeneous 
space in different and completely arbitrary coordinate systems Sa and Sb- Both of 
them apply the same definite method to transform to "optimized" systems S'^ and S'^ 
respectively. Then the systems 5"^ and S'^ thus obtained can only be related by a trans- 
formation from the eight parameter group above. 

Optimization: The metric in the optimized system should be "as close to an exact 
R-W metric as possible" within a given four dimensional volume. Since any perturbative 
treatment is performed around a R-W background, this requirement guarantees that the 
perturbation becomes as "small" as possible. One has to define the conditions of this 
optimization and to construct the "optimal" coordinate system. 

Our covariant optimization proceeds via a series of minimizations of four dimensional 
volume integrals. Expressions which vanish for the exact R-W metric are minimized by 
choosing an optimal gauge. This gauge fixing is performed as far as principally possible. 
Starting from an arbitrary system S with metric g^i,{x), one constructs an "optimal" 
system S' with the transformed metric g'^j,(a;'). At the end one can define the "optimal" 
approximating R-W metric Ij^y^x'). 

The method has the following properties and advantages. 

• It only uses the metric g^^,^, no assumptions about the origin of the energy momen- 
tum tensor are necessary. A scale factor does not show up during the procedure, 
the optimal scale factor is determined at the end. 

• Foliation of space-time is part of the procedure and is obtained in a unique way, 
again without resorting to any assumptions concerning a fiow of matter. 

• The general formulation is non perturbative. 

• The procedure can be explicitly applied in perturbation theory if the metric is a 
small deviation from an exact R-W metric. We will present the explicit formulae 
in first order of the perturbation and sketch the procedure for the second order. 

• The four dimensional volume over which the minimization is performed is arbitrary. 
For the perturbative treatment we will specialize to simple volumes. 
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• Covariance is an immediate consequence of the method. Since the variation is 
taken over all coordinate systems, it does not matter in which coordinate system 
one starts. 

• Spatial averages of arbitrary tensor fields can be naively performed when using the 
optimal coordinate system. There is no need to decompose into tensor structures 
or to restrict to averaging of scalar quantities. 

• There is no backreaction in first order. 

The paper is organized as follows. In sect. 2 we present the general non perturbative 
method. In sect. 3 this is explicitly applied in first order of the perturbation. The 
Euler- Lagrange equations for the coordinate shift which leads to the optimal system are 
formulated and solved. Boundary effects turn out to be irrelevant if the wave length of 
the perturbation is small compared to the spatial extension of the volume. In sect. 4 we 
present the rather simple extension to second order. The short sect. 5 describes averaging 
which has become trivial. Sect. 6 deals with backreaction. There is no backreaction 
in first order. In a simplified static treatment of the second order it turns out that the 
second order of the transformation is not needed, only the first order perturbation of 
the metric introduced into the second order Einstein tensor enters. We consider the 
contributions of galaxy clusters, galaxies, and stars. We give arguments why neither a 
concentration of clusters and galaxies in bubble walls which surround large voids nor 
retardation effects are relevant. We find that the ratio Pb/p between the density pi, 
mimicked by backreaction, and the averaged matter density p is small, of the order of 
10~^ to at most 10"^. To mimic dark energy from backreaction appears practically 
impossible. Sect. 7 gives a summary. 

2 General conditions for the optimal coordinates 

For an exact spatially flat {k = 0) Robertson- Walker (R-W) metric one has 

9mn = a^(^)^mn, 9m0 = 0> ^00 = "l- (2-1) 

Consider now a realistic metric with metric tensor g^^^x), and a given four dimensional 
volume. We want to define "optimal" coordinates x'^ in which the metric within this 
volume becomes, in a sense to be defined, as close as possible to an exact R-W metric 

with A; = 0. It would be impractical to combine all conditions into a single variation 
problem by minimizing an integral over the sum of appropriate squares. This would lead 
to rather complicated Euler-Lagrange equations even in a perturbative treatment. It is 
technically much simpler to proceed in steps. Each of the first four steps approximizes 
a certain property of the exact R-W metric. A scale factor does not show up in the 
conditions. The optimal scale factor a{t) associated with the given realistic metric is 
determined at the end of the procedure in step 5. 

All integrals in the four steps below are taken over the four dimensional volume under 
consideration. In step 1 the time coordinate is fixed but arbitrary, and the variation is 
taken over all primed systems which are time independent coordinate transformations 



3 



of the original one. In steps 2 - 4 we also allow time dependent transformations. These 
transformations have to respect the restrictions obtained in the previous steps. Because 
all coordinate systems are admitted in the variation procedure the covariance of the 
method is automatically guaranteed. Each step restricts the freedom of choice of coordi- 
nates more and more, at the end the coordinates are fixed as far as principally possible. 
Step 1 leads to a transversal perturbation. Usually one will start already with some 
"reasonable" coordinate system. In this case steps 2, 3, 4 become trivial. In step 5, 
finally, we define the optimal scale factor a{t) which gives the "best" approximation of 
the given metric to an exact R-W metric. Here and in the following 

if) - {fm = (2-2) 

denotes the spatial average of /(x, t). 
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Step 1 : / I ^'-"^^'''l , (fu/^^ran \ /_^,(^A^4^/ _ Minimum, 

V (4/3) ; 

where time is fixed but arbitrary, and the variation is over all time 
independent coordinate transformations. (2-3) 

Step 2 : / (^k^L=] ^-g'{x')d'x' = Minimum, 

with the variation taken over all coordinate transformations 

which respect the restrictions obtained in step 1. (2.4) 

Step 3: J [g'ooix') + if -g'{x')d^x' = Minimum, 

with the variation taken over all coordinate transformations 

which respect the restrictions obtained in steps 1,2. (2.5) 

Step 4 : j i^^g'^^{x')^ ^-g'{x')d^x' = Minimum, 

with the variation taken over all coordinate transformations 
which respect the restrictions obtained in steps 1,2,3. (2.6) 
Step 5 : Define the optimal scale factor a{t') by a^{t') = {g'J?,){t'). (2.7) 

Summation convention is always understood, also for identical lower indices. 

The meaning of the conditions should be obvious. For an exact R-W metric all the 
integrands would vanish. In steps 1 and 2 we introduced normalization factors in the 
denominator. This is necessary, because otherwise one would run into an unphysical 
minimum by a simple scaling of the metric. 

The conditions above do not fix the metric completely. They yield a whole class of 
optimal coordinate systems. This class is, by construction, independent of the system 
with which one has started. The freedom which remains are the transformations from 
the eight parameter group of rigid translations in space and time, rigid rotations in space. 
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and rigid scaling of the space coordinates. This is just the invariance group of coordinate 
transformations mentioned in the introduction which keep the form of the exact R-W 
metric intact. 



3 First order 

We consider a perturbed Robertson- Walker metric of the form 



9mn{x) = a'^{t)Smn + hmn{x), gmo{x) = hmo{x), goo{x) = -1 + hoo{x). (3.1) 

The gauge and the way of splitting into background and perturbation is completely 
arbitrary, except that the perturbation should be small. The scale factor a{t) is in 
general not identical with the optimal scale factor a{t) obtained at the end. 

Introduce the perturbed metric in the primed system into the integrands in (2.3) - 
(2.6). The primed system is connected to the old one by an infinitesimal transformation 

= 2;'^+^^. Because all the brackets vanish for the unperturbed metric it is sufficient to 
expand these up to first order in /i'^^ and ^f^. Furthermore we can put = a^{t) and 
restrict to the leading order in the denominators. Express the h'^i, in the primed system 
by the h^j^^, in the old one and the shifts In lowest order one has = 2^(t)^'", ^0 — 
— The well known transformation laws for the metric in lowest order of ^fj, are: 



= a'^{t')Srnn + hrnn + ^m,n+^n,m —'^0.{t')a{t')^oSrnn, (3.3) 

g'moi^') = 9mo(x) + Co,m + Cm,o " 2(a(i)/a(i))^„ 

= hmo + Co,m + a''it')r,o, (3.4) 

^00(3^') = i/oo(a;) + 26,0 

= -1 + /loo + 26,0 ■ (3.5) 

The transformation of the averages {g'^^){t'), which now refer to a different time t', is 
most easily obtained by writing d^y' — 5{y'^ — t')d^y' . This results in 

{g'J^) - S^(^) + (Y + ^e.^)+2a5(eo-(eo)) (3.6) 
= 5^(0 + + 255(6). (3.7) 

In this way we obtain the integrals which have to be minimized. 
Step 1: 

We have to minimize 
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/( ^'~'^|-g;^'>'"'" )'v^^^'= (3.8) 

f 1 f 1 2 ^ ^ ^4 



with respect to while keeping arbitrary but fixed. Here and in the following it is 
irrelevant whether we consider the expressions in the old or in the new system. In all 
brackets the leading terms cancel, changes in the boundary of the volume only contribute 
to higher order. 

From a variation 6^rn in the interior we obtain the Euler-Lagrange (E-L) equations 
We use the usual decompositions for /i^„ and hmo' 

^mn ^ (t){^^mn ~l~ .^,mn ~l~ C^^^ -\- Cn,m ~l~ ^mri\i (3.10) 

hmo = d{t)[F,m+Gml (3.11) 

with 

Gm,m — 0, Drnn — -Dnrm 0, Drnn,n 0, Gm,m 0. (3.12) 

A special solution of the E-L equations (3.9) is then (here A = dndn) 

e^)(x, t) = -a'it){U-\A - 2^^^o),m +lB,m (x, t) + C^(x, t)}. (3.13) 
Introducing this into the transformation gives h'mn as in (3.10), where now 



A' ^ A- 2^fo, B' = -A-^(A - 2^eo), C'm = 0, D'^n = Dmn. i.e. 

a[t) a[t) 

h'mn = a^it')[idmdn-5mn^)B' + Dmn]- (3.14) 

We have A' + AB' = 0, which implies that h'mn is transversal, h'mnm = 0. 

The solution ^j^^ in (3.13) is not unique, because neither the operator nor the 
decomposition (3.10) is unique. To see this more explicitly, let c be a constant and </?, 
functions with A(p — At/j — 0. Then one can replace A — >■ A+c, B — >■ B—cx'^/2—'^, Cm — >■ 
Cm - f,m,Dmn ^ ^mn + '^'^mn +}P,mn, wlthout changing hmn- In particular one can 
always remove a constant D„in = D,nn by putting it into B,mn with B = DijX'x^ /2. We 
assume that this has been performed if necessary. For a detailed discussion we write the 
most general solution of (3.9) as 

e„.(x,t)=e^)(x,i)+77^(x,t). (3.15) 
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Introducing this into (3.8) leads to a variation problem for r]m'- 

J {Vm,n +Vn,m -^{riin )^mn + ^ mn " g (^'n)(^mn) = Minimum. (3.16) 

A variation in the interior results in the homogeneous equations associated with (3.9), 
i.e. 

Vmmn~^Vnnnn 0- (^■■^'^) 

Two tasks have to be done. Firstly one has to consider variations which also involve 
changes at the boundary. This will yield a special solution for ?7m(x, t) which fulfills the 
boundary conditions. Secondly one has to classify the whole set of solutions, in order to 
determine the remaining freedom. 

Concerning the first point it will turn out that the solution is concentrated in a strip 
along the boundary, with an extension of the wave length of the perturbation h'mn- If the 
wave length is small compared to the extension of the volume it is therefore irrelevant. 

As for the second point we will find that the freedom consists in rigid translations, 
rotations, and scaling, at this stage still with an arbitrary time dependence. 

If these results appear obvious, one may skip the following lengthy derivations and 
proceed directly to (3.24), (3.25) at the end of step 1. 

We now treat general variations rjjn ^ f]m + ^Vm which also involve changes at the 
boundary. This poses a delicate problem, because rjm has to fulfill the homogeneous 
E-L equations in the interior. Therefore it would not help to consider the boundary 
terms from the partial integration because rjm cannot be chosen completely free at the 
boundary. We therefore proceed in the following way. 

We know that 1]^ has to fulfill the homogeneous E-L equations. Consider therefore 
a complete set of solutions r^j^l of (3.17), expand rjm = J2aCaVm^ introduce into the 
variation problem (3.16). Differentiation with respect to the coefficients leads to the 
linear system of equations 

+ / + V^nin ' ^ (r?H)^mn) (/^'mn " l{h' j,)S^n)d'x = 0, with 
^'"^^ = / ivSn + - l{vS)SmnWSn + vl?}m ' l{vfj)Smn)d'x. (3.18) 

We may assume that I^"^^ is diagonalized, and normalized such that it has eigenvalues 
and L"^ only, where L is some length introduced for dimensional reasons. Eigenvalues 
belong to solutions of (3.17) which in addition fulfill 

vSn + V^:in-l{vS)^mn^O. (3.19) 

Contributions of this type solve (3.18) trivially for arbitrary Cq, and may always be added. 
The only solutions of these equations are rigid translations, rotations, and scaling, at 
this stage still with an arbitrary time dependence. To show this formally, we first observe 
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that (3.19) implies ry^\nk + nn,mk — (^^ ^^^^ point it becomes apparent why we used 
the average {g'u/^) instead of g[^/?> in (2.3), (3.8)). Applying (3.19) again in order to 
exchange m, k and n, k, respectively, implies —'2.7ikinn — 0> i-®- ^ second derivatives 
vanish. Thus r^j^^ can only contain constant and linear contributions. Introducing a last 
time into (3.19), one finds that the linear part is restricted to a scaling and a rotation. 

The coefficients which refer to the non trivial solutions associated with the eigenvalues 

are uniquely fixed, namely 



~h I ~ \{h'n)5^n)nndA, (3.20) 



with rin the normal vector at the boundary. The term ~ {Vi,i )^mn in the first line does 
not contribute, in the second line we performed a partial integration, making use of 
h' =0 

We give examples with plane waves in 2;— direction. Consider first a gravitational wave 
h'mn = Knn COS kz, with Jin = — /i22 = h = coust. as the only non vanishing components. 
This implies h'jj = 0. It is convenient to choose the volume as a cylinder with radius po 
and < z < L. Introduce cylindrical coordinates p, </?, z, together with the corresponding 
unit vectors e^^ = (cos<^, sin(y9, 0), e^^^ = (— sin (/?, cos (/?, 0), — (0,0,1). Because of 
h'mz = the only surface which contributes in (3.20) is p = po, where Un = el^\ One 
has el^hmn^n^ = hcos2(p, e^hmn^^'' — — /isin2<^, e^hmn&n^ = 0. This implies that 
one only needs to consider solutions r/^ with a corresponding structure, such that the 
surface integral in (3.20) is non vanishing. Define the vectors in cylindrical coordinates 
"H^ — Vm^m^, 1 — PiVi -2- Then a basis of relevant solutions of the free equation (3.17) is 
obtained by an ansatz of the form 



/ f (p) cos 2ip cos kz \ 

g{p) sm2Lp cos kz , (3-21) 
\ /i(p) cos 2(^ sin A;^ / 

where the index (7) denotes the components in p, (p, z-direction. Obviously there is no 
need to consider any other contributions of cosn^p, sinnp. Either they cannot fulfill 
(3.17), or they give a vanishing Cq. If L is a multiple of the wave length A = 27r//c, other 
modes in kz will also not contribute to c^. 

Introducing into (3.17) leads to three coupled differential equations for the compo- 
nents p, (p, z: 



2f" + 2^ - 6^ - k^f + 2^-6^ + kh' 

P P P P 



9" 



9' 



94 - eg - 2^^ 
P^ P 



4 - 2k'l 

p2 



h" + 



h' 



P 

9i 



4\-2k^h-k[f+^-r^- 
P P P 



0, 
0, 
0. 



(3.22) 
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A (non orthonormalized) basis for the regular solutions is 




V 






(3.23) 



with Im = Imikp) the modified Bessel functions with argument kp. 

There is no need to go into more details, the qualitative behavior can be read off 
immediately. The solutions r]^^ will be combinations of modified Bessel functions Im{kp). 
These are monotonically increasing and behave fike {kp/2)"^/m\ for small kp, and fike 
exp(fcp)/\/27rfcp for large kp. If the radius of the averaging volume is large compared to 
the wave length, i.e. if kpo ^ 1, the modified Bessel functions as well as the shifts i]^"'^ 
are only relevant in a small strip along the boundary p = pq, with an extension of the 
order oi X — 27r/k. 

One can normalize /'"^l by multiplying Tyj^l in (3.23) by factors ~ Lexp{—kpo). 
For the orthonormalized solutions one gets Car^j^l ~ (h/k)^ po/ pex.p{—kpo) exp{kp), i.e. 
CaVm^ ^h/kaX the boundary, independent of po and L. Therefore, if desired, one may 
perform the limit po ~^ oo, and/or L — )■ oo and one is sure that the solution stays finite. 

If the boundaries for z are less convenient, i.e. if L is not a multiple of A, the situation 
is slightly more complicated. Instead of (3.21) one needs a superposition of terms with 
cos iTmzjL and sin 2T:nz/L, n = 0, 1, 2, • • •. The coefficients in front are only sufficiently 
large for n with n ~ A;L/27r, the coefficients for smaller n decrease in the same way as the 
extension of the modified Bessel functions increases, again the solution is concentrated 
at the boundary p = po- 

In our next example we consider a wave with B' = B cos kz, which implies h'mn ~ 
hjnn cos kz, with hii — — h — Bk^ as the only non vanishing components. Again, 
for a wave, (h'jj) ^ 0. We now have e(^)/i„„e(f) = h, e(^)^™„e(f) = e(^)/i^„e(f) = 0. 
Therefore one can use an ansatz like (3.21), with cos2(/9 and sm2(p replaced by 1. The 
further calculation proceeds as before with an analogous result. 

The previous findings are very convenient. For a large volume one can neglect the 
boundary effects, or alternatively, apply the result of the transformation only within a 
shghtly smaller volume. 

After step 1 we have transformed to a system where 

^« = a" {t') 5mn + h'^^ = a' {t') [5mn + {dmd„ - Srnn^) B' + D'^^] , 

i.e. A' + AB' = 0, = 0, ^ h'^,,^ = 0. (3.24) 

The freedom of remaining transformations in the new system, due to the freedom in the 
homogeneous solution 77"*, is now restricted to 

r = U/a\t)=b"'{t) + S{t)x"' + [iu{t)xxr, 
^0 completely arbitrary. (3.25) 
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The functions b"^{t), S{t),u!''{t) still have an arbitrary time dependence at this stage. 

Before proceeding we will assume that the above transformations have been per- 
formed, such that (3.24), (3.25) hold. These conditions have to be maintained in the 
following steps. We now drop the primes for the optimized system obtained after step 
1, and use the prime for the optimized system of step 2. 

Step 2: 

According to (2.4) we have to minimize 



/( 



V(-^'00)(^'n/3). 



^-g'{x')d^x' = J a(t)(/i„o + eo,m+S'r,o)Vx, (3.26) 



where is now restricted to the special form (3.25). The functions in the minimization 

procedure are b"^ (t) , S (t) , u'^ (t) , and ,^o(x, t). 

Introducing the decomposition (3.11) and writing the terms with 6™(t) and S{t)x"^ 
in as gradients one obtains 

h'mo = hmo + ^o,m+a\t)C,o (3.27) 
= {Co + a{t)F + U\t)p{t)x^ + S{t)x^/2]},^ +a{t)G^ + a\t)[uj{t) x xf. 

Obviously one cannot determine ^o, S{t) separately because they only enter in the 
combination 

^o{^,t) = ^o{x,t) + a{t)F + a\t)[b\t)x'' + S{t)xy2]. (3.28) 

Therefore we may, at this stage, choose b''{t) and S{t) arbitrary and vary only ^o(x, i) 

and u'^{t). 

Only for simplicity and the sake of obtaining more transparent formulae, from now 
on we specialize to the case that the volume is a sphere of radius ro = ro(t) around 
the origin. Obviously this condition is coordinate dependent, nevertheless it does not 
destroy the covariance of the procedure. In any other system, obtained by an infinitesimal 
transformation, the volume is a distorted sphere with an infinitesimal modification of 
the boundary. This only introduces irrelevant boundary effects of higher order. 

The E-L equation and boundary condition {rim denotes the normal vector at the 
boundary of the three dimensional sphere) for become 



A|o = 0, (3.29) 
(^O)m +S(t)Gm)^m — at the boundary. 

We dropped the boundary term a'^{t)[uj(t) x x]"^nm which vanishes for a sphere. 

The situation is similar to that in the first step, but simpler because we now can 
use the Neumann type boundary conditions for ^q. As an example consider a plane 
transversal wave = GSmx cos kz. This gives Gmn-m — GsinO cos (p cos{kro cos 9) . The 
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solution of (3.29) can then be expanded as .^o = G J2^iCi{r /toYP/- {cos 6) cos Lp, with 
ci = -[a{t)ro{2l+l)/2l\l + l)] £ sin^ QPlicosQ) cos{krocosQ)dQ. likro = 2nro/X > 1, 
the factor cos(A;roCOS©) oscillates rapidly, therefore the integral is only important if 
this oscillation matches with the oscillation of the Legendre polynomial, i.e. if / ~ 
2A:ro/7r. This is large, therefore the factor (r/rg)' is only relevant near the boundary. 
In the interior the solution is essentially zero. We may therefore replace (3.29) by the 
corresponding Neumann problem with = which has the unique solution = ~'?"(^)) 
where r is constant in space but may depend on time. Therefore (3.28) gives 

eo(x, t) = -T{t) - a{t)F - a\t)[h\t)x^ + S{t)x' /2]. (3.30) 

Next we vary ijj^{t). Because this is a function of t only, the x-integration remains and 
we obtain 

We used some simplifications for the case of a sphere, and the result just obtained for 

lo, 

After applying the transformation in (3.27) with the results (3.30) and (3.31) we have 
transformed to a system where 

h'^^ = a{t')G'^,Wii\iG'^ = Gra + Mi'){^'^Gn-x"G^)x''/2, 
i.e. F' = 0, ^/i'^o,m=0. (3.32) 

The remaining freedom for coordinate transformations in the new system is now 

r = U/a\t)^h^{t) + S{t)x^ + [ujxxY', 
^0 = -^^^T{t)+o?{t)ih\t)x'' + S{t)x^/2]. (3.33) 

We are no longer free to perform different rotations at different times, — const., while 
translations h^{t), scalings S{t), as well as r(t), can still be time dependent. 

Again we assume that the transformations have been performed and that the follow- 
ing steps respect (3.32), (3.33). The primes will be dropped. 

Step 3: 

We have to minimize 

j (^'oo(^') + lf^J-g\x')d''x' = j d\t)(hoo + 2Uofd''x. (3.34) 
Using the result (3.33) for ^o, the bracket becomes 

/ioo + 2^0,0 = /ioo - 2^{T(t) + U\t)ih\t)x'' + S{t)x^/2]}. (3.35) 

Variation of T{t),b''{t), S{t) gives the equations (3.36) - (3.38) below, in which we al- 
ready performed the spatial integrations where possible. In the derivation we used some 
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simplifications which hold for a sphere, i.e. (x™) = 0, {x^x^) = 1^x^)5^^ /?>. For a sphere 
of radius ro one has (x^) = 3ro/5, {{x^Y) = 3ro/7, and, in the denominators below, 
{{x'^f) - {xY = 124/175. 

{h,o)-2f{t)-{x')^{a'{t)S{t)) = 0, (3.36) 

{x^hoo)-l{x')^^{a'{t)b^it)) = 0, (3.37) 

{x'hoo)-2{x')f{t)-{{x'f)^{a\t)S{t)) = 0. (3.38) 

The solutions are 



^ ' " 2 - ' 

|(a^W*"W) = (3.40) 

The function r(t) is now fixed up to an additive constant r, while b'^{t) and S{t) are 
fixed up to two integration constants 6"*, and 5", E, respectively. 

b-^{t) = b-{t) + b-^ + rj S{t) = S{t) + S + Ej (3.42) 

with b"^{t) and S{t) obtained from integrating (3.40) and (3.41). 

After having performed the transformations of step 3 one has the equations {h^o) = 
(a;"^/iQo) — i^^hQo) — ■ The remaining freedom for transformations is now 

r = ura\t) = b^+sx"' + [uxxr+{r + ^x^) j 

^0 = _^Q^^ + /3V + ExV2, (3.43) 

where b"\ S, ou'', r, S are all constant. As before we assume that the transformations 
have been performed and that the form of further transformations is restricted to (3.43). 

Step 4: 

Before proceeding with this step we look at the meaning of the constants /3™ and S. 
Obviously /3™ describes an infinitesimal boost (slightly modified because a^(t) 7^ const.), 
while E describes a time dependent scaling. A transformation with (3.43) gives a space 
dependent contribution 255 [P'^x'^ + Ex^/2] 6jnn in g'mm which illustrates that boosts do 
not leave the standard form of the R-W metric invariant. This contribution is, however. 
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suppressed by the factor a, therefore a minimization along the previous hues appears 
ineffective. We therefore choose condition (2.6) in step 4, and minimize 

= / 53(t)(/ioo,o+~^(/3' + Sa;'=)/ioo,fc)V 



We are left with an ordinary minimization problem for the four constants P'' and S. 
Define the integrals 



Irnn 



J 
J n 




hooTmhoo:0 d^x. (3.45) 



The linear system for (3'' and E then becomes 

7E + 7„^" = -J 

Im^ + Imn^^ = -Jm- (3.46) 

This can be easily solved after Hqq is specified. Practically it is even simpler. For a large 
averaging sphere the volume integrals are rotation invariant, which implies 1^ = Jm — 
and Imn = Smnlii/^, leading to ^"^ = 0, E = - J/7. 

In the transformed system we will then have J' = 0, J'„^ = 0. 

The nature of step 4 is somewhat different from the steps before. While steps 1, 2, 3 
only require that h^i, is small, step 4 requires in addition that time variations are small 
compared to spatial variations. 

Our optimization procedure has now come to an end. The integration constants 
E are fixed. After step 4 the only allowed transformations which remain are those 
which keep the form of the R-W metric invariant. 



r = U/d''{t)^b"' + Sx"' + [ujxx]'^, 
e = -eo = T, (3.47) 

with h^,S,u^,T all constant. 
Step 5: 

In the last step we determine the optimal scale factor a{t). In (3.24) we obtained the 
form g'^^ = a^(t)[(l — AB)dmn + -B,mn +-Dmn], where we denote the first order result for 
the optimized metric by an index (1). According to (2.7) we define 

a'{t) ^ {g^/m) = - \{^B)]. (3.48) 
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Eliminating a'^{t) gives g^^l = a'^{t)[{l - AB + |(A5))5^„ + B,^n+Dmn], which can 
finally be written in the form 

g^Ji = a\t)[5,r^n + {dmdn - Smn^)B('^ + (3.49) 
with = S - (AS)x76, D« = Dmn- One now has {h^^J = -2a'^{AB^^^) = 0. 

The whole procedure is simpler than it appears. Usually one starts already with an 
ansatz for the fluctuations in a "reasonable" coordinate system. Then steps 2-5 may 
become trivial from simple symmetry arguments, i.e. one is already using the optimal 
coordinates and no further transformations are necessary. 

Summary of the first order transformation 

We found that the transformed metric g'^J after the flrst order optimization has the 
form 



9ml = a\t)6mn + hmn = a\t)[6mn + {dmdn - SmnA)B + Dmn], (3.50) 

9mo = ^mo = a{t)Gm, (3.51) 
gS^ = -l + Zioo- (3.52) 

The approximating R-W metric g^^, is obtained by dropping the perturbations B, 
G„i, /loo- This provides a natural basis for splitting into background and perturbation, 
and for performing perturbative calculations. 

The quantities which could not be removed by gauge transformations nevertheless 
share some properties of the unperturbed metric. Several relations can be derived by 
appropriate partial integrations and some simplifications which hold for a sphere: 



hmn.n = 0, {K) = 0, = 0, ((x™x" + xH^''/2)hmn) = 0, (3.53) 

hmOrm = 0, {x^'hmo) = 0, {XmhnO — Xnhmo) = 0, (3.54) 

(/loo) = 0, (x"^/ioo) = 0, (a;2/ioo) = 0, (3.55) 

j a^it) I ^1 I ^oo,m /ioo,o d'^x = 0. (3.56) 

The relations (hmm) — 0, (hoo) — 0, will lead to the absence of backreaction in first 
order. 



4 Second order 

We will see in sect. 5 that the second order of the transformation to the optimal gauge is 
not needed if one neglects time derivatives in the background metric and in the perturbed 
density. Nevertheless, from a principle point of view, it is instructive to show how the 
procedure can be extended to second order in a quite simple way. 
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Rather than trying an approach by brute force, one should proceed in an iterative 
way. Determine the shift which leads to the optimal system in first order. If one uses, 
instead of (3.2) - (3.5), the exact transformation formulae, or at least, considers up 
to second order, one obtains a perturbation of the form h^i, — h^^} + h^j, in the optimal 
system. All observers, irrespective of their original gauge, have ended up with the same 
/i^^j, up to the remaining 8-parameter group of rigid transformations. For the higher 

order contributions hj^i, which have not been optimized this is not the case. Unlike h^^^ 
they do depend on the original gauge. An explicit calculation would be tedious. We will, 
however, see that h^i, does not contribute to backreaction, therefore there is no need to 
calculate it. 

We now start the second order calculation with the metric h^^ = + f^fii^- Within 
the brackets in (2.3) - (2.6) we have to expand the metric gf'^^^ up to second order in the 
perturbation and to perform the minimizations in the four steps. There is a considerable 
simplification due to the iterative procedure. Because /i^^j is the solution of the first 
order problem there are no first order terms in Consequently needs only to be 
considered in lowest order, i.e. one can again apply the simple transformation formulae 
(3.2) - (3.5). 

It is convenient to include the factor \/ —g' into the squares by writing 

1,}^ 



^/^ = a^[l + ^(^-/^S)] . (4.1) 

The scale factor a has been replaced by the optimized scale factor a of the first order. 
The brackets have now to be considered in second order of the perturbation. Due to the 
properties (/i^?^'') = and (/iM ) = 0, there are no corrections from the averages in the 
denominators of steps 1,2. 

The shift is a superposition of two terms, = + C^*^"*- The first term, 
corresponds to the solution of the E-L equations with h^.^ replaced by hy,y. This part 
transforms hp^y into w-"^ in the same way as previously it transformed /i^^ into h!-^. 

Therefore h'^'^l shares the properties (3.53), (3.55), i.e. (Ji^^^) = 0, (h^QQ) = 0, consequently 
it will not lead to any backreaction, and one does not need ^j^^ and h^^} . 

The second contribution, is due to the quadratic terms in h^^} in the brackets 
in (2.3) - (2.6), which now have to be inserted into the E-L equations. This part could, 
in principal, lead to a backreaction. As mentioned before, we will however see that this 
does not happen if we drop the time dependence. 



5 Averaging 

In first order of the perturbation we finally obtained an optimal coordinate system with 
a metric of the special form (3.50) - (3.52), and the approximating R-W metric 
obtained by dropping the perturbations. One can now simply define spatial averages of 
arbitrary tensors in the naive way as in (2.2), i.e. 
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The freedom in the choice of coordinates is restricted to the eight parameter group of 
global transformations described by the parameters bm, S,Uk,T, and these transforma- 
tions commute with the operation (5.1) of averaging. 

It is not necessary to restrict to averaging of scalar quantities or to decompose into 
invariants (usually only with respect to purely spatial transformations) and perform the 
averages for the latter.. In fact such an approach can be problematic. The expansion 
tensor Qmn = gmn,o /2 is often decomposed into invariants, and quantities which enter 
linearly and quadratically are averaged separately (see e.g [2]). This procedure can 
lead to quite strange consequences like negative averages of positive definite expressions. 
Nothing of this kind can happen in our case. 

For the special case of the metric we note that (Qfiu) is not necessarily identical with 
g^^. Due to (3.53), (3.55) the relation is, however, true for qqq and for gmm in first order, 

i.e. ((^oo ) = ^00 ^nd ((7^^) = gmm- These are the relevant quantities for backreaction 
which we discuss now. 



6 Backreaction 

In the same way as the perturbed metric, g^u = g^^ + Sg^^,, we split the Einstein tensor, 
G,_iu = G^y + 5G^y, with G^y the Einstein tensor associated with gfj_y. Using the Einstein 
equations (we include a cosmological constant A) for G^^ one has 

G^y = (G^y) = {G^y) - {6G^y) = k{T^,) - A{g^,) - {SG^,). (6.1) 

The first two terms, k{T^u) —-A{giJ.u), describe the equations which one would expect from 
the averaged energy-momentum tensor and metric, the third one, —{6Gf^y), is the devi- 
ation from this, i.e. the backreaction. The essential quantities associated with density 
and pressure mimicked by backreaction are Kph = —{SGqo) and upb = —{SG^)/3. In first 
order of the perturbation one has (indices are raised and lowered with the background 
metric g ) 



^Cffi = i(/^^^.,-/^i,^)-2^/^^ + -/^i,o-2^/^^„, (6.2) 

^ Qi Qi Qi 

Q Q Q 

• 2 •■ 

-3(^ -f 2-)/zoo - /ioo,i -3-/100,0 • (6.3) 

If one inserts h^^} and uses the properties (3.53) - (3.55) one observes that there are three 
^(1) 



types of terms in 6Gqq and 6G^"^: 
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a) terms which vanish, 

b) terms where the spatial average vanishes, 

c) terms which are spatial derivatives and can be written as surface contributions in 
the integral which are irrelevant for large volumes. 

Therefore {6G'£) = {SGg^"") = 0, there is no backreaction in first order. This is in 
fact a property which one expects from any reasonable lowest order averaging prescrip- 
tion, where positive and negative contributions of the fluctuations should cancel. 

We now come to the second order. Here one has two types of contributions. The first 
one arises from introducing the second order correction of the metric into the first 
order correction {5G^^^) of the Einstein tensor. In general this does not vanish because 
does not fulfill the properties (3.53) - (3.55). For an estimate one can, however, 
use a simple static approximation (we will comment on retardation below) which is 
appropriate for the present day universe. The peculiar velocities of galaxies are small, 
with v/c of the order of 10~^. We also do not have sizeable perturbations with extremely 
short wave length which could contribute large derivatives. Let the extension in time 
of the averaging volume be small compared to the Hubble time. Then the variation of 
h^i, in time, caused by the slow motion of matter, as well as the weak time dependence 
of a{t), can be neglected in comparison with the variation in space. It is convenient 
to fix a{to) = 1 for the present time to- The surviving parts of {SGqq) and (^G^^™") 
in (6.2), (6.3) only contain spatial derivatives, i.e. the averages vanish up to irrelevant 
boundary terms. This holds irrespective of the special form of and is, of course, very 
convenient. There is no need to calculate h^^^. 

The second contribution arises from introducing the first order correction h^^^J of the 
metric into the second order correction {5G^^^) of the Einstein tensor. The rather lengthy 
expression for 5G^^^ can be found in Wetterich [TU] in the approximation that derivatives 
acting on the background metric are neglected. It will not be written down here. We 
only give the two quantities needed for the effective density and pressure, anticipating 
hmo = 0, neglecting the time dependence in the perturbation, and performing some 
spatial partial integrations. 

= {^hooK;i-lhooh^^,,,+^K,^,t+lh^=h,,,l-^h^^^ (6.4) 

(<5G(^)-) = {-hooho,l +lhooh''n, +1^1^,1 +\h''h],,, |6.5) 

In order to estimate h^^, = h^^} from the sources we use again the simple approximation 
above, i.e. neglect any time dependence. This is essentially the model of Wetterich [TO] . 
It has the advantage that it is transparent and leads to explicit formulae. 

Consider a dust universe with p only weakly time dependent and p = 0. The solution 
for the perturbed metric is only needed in lowest order and most conveniently first derived 
in the harmonic gauge (marked by a hat) 

K,.= h:„, (6.6) 
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where the perturbed Friedmann equations have the simple form 

h,,,P=-2Ki5f,,-^5f^g^,). (6.7) 

In our simple model 6Tqq = 6p is time independent while all the other components of 
ST^^ vanish. This implies 

^^00 = — /t5p , furthermore = (J^/ioo, , ^mo = 0. (6.8) 

One can transform back from the harmonic gauge to our optimal (transversal) gauge by 
a purely spatial shift ^"^ = 

^00 = ^00) ^" 

Introducing (6.9) into (6.4) 
and one obtains 

Pb = 
Pb = 

The pressure term coincides with the result of [10], the density term has a factor 7/4 
instead of 9/4 in [10]. This shows that the correction is smaller in our "optimal" gauge 
as one would expect. It also illustrates that differences between reasonable gauges are 
small. 

We now consider a more specific simple model. We start with a hierarchy of clusters, 
composed of galaxies, composed of dark matter and stars. Subsequently we will also 
discuss the modifications due to the presence of voids. 

Consider first a space filled with homogeneous spherical objects of radius L, density 
p, and mass m = AirpL^/S, which are roughly uniformly distributed in space at positions 
Tj. The average density is denoted by p. The corresponding density fluctuation is 

Sp{r) ^ p(r) - p = 0(L - |r - r,|) - p, (6.12) 

i 

such that the average Sp vanishes. To find a useful approximation to the corresponding 
/ioo(i") we first define a distance D by the requirement that the average of 6p vanishes 
within a sphere of radius D around a source, i.e. D^/L^ = p/p. This D is roughly half 
of the average distance between the spheres. For well separated sources one has D ^ L. 
A solution of (6.8) within this sphere (chosen around the origin for simplicity) is then 

r 2 2 r 3 2 

hooir) = k{p (— - j)e{L - r) + p — e(r -L)+pj + pc}. (6.13) 

This is just the well known potential of a uniformly charged sphere in a constant back- 
ground. The constant c has been introduced in order to achieve /^<^ hood^x = 0. It is 
of the order L^/D and will turn out to be irrelevant. Of course (6.13) is not an exact 



—A ^hQQ,jn /2 which leads to 

, = (C - A-'d^d^)hoo, h^o = 0, 5T^, = 5%,. (6.9) 

, (6.5), all terms can be expressed by (/looA/ioo) = —i^ihoo^p), 

--{SG^S) = lihooSp) (6.10) 

K 4 

-^{SG^^^n = -^{hooSp) = -^P.. (6.11) 

6k VI IV 
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solution of (6.8), (6.12). The spheres of radius D around the sources overlap in some ar- 
eas and leave empty regions elsewhere. Furthermore /iqo ~ KpL"^ !?> 7^ at the boundary 
of the circle, thus introducing boundary contributions there. The average of the ficti- 
tious mass distribution associated with these corrections vanishes, and the location is a 
distance ^ D ^ L away from the sources where /iqo enters in (6.10), (6.11). Therefore 
the resulting corrections to /iqo are suppressed compared to (6.13) and can be ignored. 
The average (/ioo) will already vanish when taken over regions involving only a modest 
number of sources. This implies that also (x™'/ioo) = and (x^/iqo) = 0, as required in 
(3.55). 

We are interested in the ratio pb/p = (7/4) (/iqo ^p)/p- When calculating {hoo6p) one 
can drop the constant —p in (6.12) because (^00) = 0, i.e. one can replace Sp{r) by p(r). 
Only the regions |r — r^l < L where p(r) 7^ contribute in the product. This results in 
a considerable simplification. It is now straightforward to calculate (/ioo<^p) from (6.12) 
and (6.13). Dropping corrections which are suppressed by higher powers oi L/D one 
finds the following result which can be written in various useful ways. 

Pb/p = -{hooOp}/p = —npL = Yo'^P = ^i^j; = Y^^ooiL). (6.14) 



Before proceeding we give a more careful justification for our neglect of time dependence. 
It has been argued by Kolb, Matarrese, Notari, Riotto [11], as well as by Bochner [T2] . 
that retardation effects, though irrelevant for nearby sources, may become important 
when taking into account contributions from distant regions. We thus look at retardation 
effects. Let D be the average distance between the sources, v c their average velocity, 
and define a distance D such that {v/c)D <^ D <^ D. Retardation is negligible as 
long as the sources can only move a small fraction of D within the retardation time, 
i.e. for {v/c)\r — r'| < D -C D. In the additional contribution hoo, where retardation 
might become relevant, one has the retarded solution hQo{r,t) = [k/Att) J 6p{r + r",t — 
r" /c)d^x" /r" , where the integration is only over the region r" = |r — r'| > {c/v)D ^ D. 
Because r" is large compared to the distance of the sources, and 6p homogeneous on 
average, this expression is independent of r, i.e. it would contribute a non trivial, time 
dependent but space independent h^Q. But this is impossible in our gauge because it 
would contradict the condition (/iqo) = 0. 

Here it becomes clear that the question of backreaction is not only a problem of using 
a reliable approximation for the distribution of matter and a consistent mathematical 
treatment. It is also crucial to connect observations with statements about the metric, 
in particular about the scale factor which describes the expansion. If observations refer 
to something like our "optimal gauge" there is no back reaction from retardation. If, on 
the other hand, they refer to some "bad" gauge in the past, there might be backreaction 
effects. Observations use light which essentially moves through a space which expands 
according to the average density p, the corresponding average fluctuation 6p vanishes. 
Therefore we don't expect a sizable additional contribution hoQ. 

If the universe would be only built up from clusters, and if clusters would be homoge- 
neous objects, (6.14) would be the final result and one should use it with the assignment 
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P ^ PCi L ^ Lc, D — )• Dci rn — )■ nic- We keep this contribution and next consider the 
additional effect that clusters are made up from galaxies. Because the average effect of 
clusters has already been considered, the additional density fluctuation within a cluster 
is described by an expression like (6.12), where now rj denotes the position of galaxies, p 
has to be taken as the density within a galaxy, and p is replaced by the average density 
Pc of a cluster. The product of /iqo and 5p in {hoo6p) does not contain mixed terms be- 
tween the expressions for clusters and for galaxies. The reason is that the terms referring 
to galaxies contain contributions which are located at distances ~ Dq which are small 
compared to the extension Lc of the clusters. Therefore one may replace them by their 
averages over the cluster which vanishes. Within one cluster we thus obtain again the 
contribution (6.14) where now the parameters are those for galaxies, and the average as 
well as the average density refer to a single cluster. Because of p = pcVc/Vtot one has 
(/loo ^p)c/Pc = {hoo ^p)/p- If only a fraction i] of galaxies is located in clusters and the 
fraction 1 — 77 roughly uniformly distributed, the expression for the clusters is multiplied 
by a factor 77, while rj cancels in the two contributions for galaxies. 

Finally consider that galaxies are made up of stars and dark matter. Because only 
a fraction of matter is made up of stars, while dark matter is assumed uniformly dis- 
tributed, there is a factor Qb/Qm- Thus the contributions of clusters, galaxies, and stars 
add up to 



Wb = Pb/pb = - — - (6.16) 

This expression has the same structure as the result in [10], but with the definite factor 
21/407r = 0.17 in front, which, of course, should only be considered as a rough estimate. 
Our ratio Wf = pj/pf = —1/21 is slightly different from the ratio —1/27 obtained in 
[To] , which is due to the different gauges. 

Galaxy clusters and galaxies are not distributed uniformly, but there exist large 
voids, surrounded by bubble walls. Because this fact plays some role in the discussion 
on backreaction, let us discuss the implications. Again we use a simple model which 
shows the essential features and can be treated explicitly. We describe the bubble wall 
by a uniform distribution of matter in a shell of radius R, thickness 2L, density p, and 
total mass m. For L <^ R, which we assume, one has p = m/SiTR'^L. Let the bubbles lie 
close together with only little space left in between. The density fluctuation for a bubble 
centered around the origin is then 

Sp{r)=pe{L-\r-R\)-p, (6.17) 
with p = p 6L/R. The solution of (6.8) with (/iqo) = becomes 



ioo(r) = Kp\ - IrLQ{R -L-r) 



+ ,'^!±i!_!fli_:!_(^)e(L-i-fl|) 

^ 2 5 6 3r ^ ^ ' 
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+ {^-^I^±^^-^^RL)Q{t-R-L) + ^t' + cV (6.18) 

V O -TL 

with c = 0{L'^). In contrast to the case of matter located inside a sphere, one can here 
perform the hmit L — > 0, while keeping the total mass fixed. This results in 

Following the same steps as before one arrives at the back reaction due to bubble 
walls, 



{pb/p)BW = -^{hoo Sp)/p = —upRL = —up R^ = ^/^^ = -^hooiR). (6.20) 

We next estimate the various contributions by inserting some standard values. 

For clusters we use rric = lO^^ruQ = 2 ■ W^^g, and Lc — 5Mpc — 1.5 • lO^^cm. 

For galaxies we insert the values for the milky way, rriG = lO^'^rriQ = 2 ■ lO'^^g, and 
Lg = lOOkpc = 3 ■ lO^^cm, where both numbers include, of course, both baryonic and 
dark matter. For the contribution of the stars, finally, we insert the values of the sun, 
m5 = mo = 2 • lO^^g, L5 = = 7 • 10^°cm. 
For the ratios mass over radius one obtains 



mc/Lc = 1.3- 10'^^ gcm~^ (6.21) 
ma/LG = 0.7 ■W'^^gcm-^ (6.22) 
ms/Ls = 3- 10^2^ cm-\ (6.23) 

We further use k, = 2 ■ 10^^"^ g^^cm and Qb/Qm = 0.17. For all three cases one has 
|/ioo('")| < /ioo(O) < 0.3 • 10""^, furthermore /iqo is smoother than 6p. Serious doubts 
concerning the validity of perturbation theory appear inappropriate. Nevertheless one 
can often read the argument that a perturbative expansion would be inadmissible because 
Sp/j) is large. One should take a closer look at this "argument". In our model Spmax = 
Sm/AnL^, p ^ Sm/AnD^, therefore indeed Spmax/p ~ D^/L^ is large. But there is no 
reason to panic, because this quantity has nothing to do with the perturbation expansion. 
The relevant dimensionless quantity which enters is not Sp/p, but nm/L which is small, 
less than 3 • 10~^ in all cases! 

The three contributions in (6.15) become 



{Pb/p)c = 77 -0.4 -10-^ (6.24) 
{p,/p)g = 2-10-6 (6.25) 
{pb/p)s = 1.5 -10-^ (6.26) 

All these corrections due to backreaction are very small. 

Let us look at some other stellar objects with a larger ratio m/L which might be 
relevant. 
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For white dwarfs the ratio m/L is about a factor of 30 -j- 100 larger than that of the 
sun, the number of white dwarfs is estimated as about 10% of all stars. Therefore their 
contribution is also unimportant. 

For neutron stars and black holes one has m/L = 6 ■ 10^^(7 cm~^, which is a factor of 
2 ■ 10^ larger than the ratio for the sun. If rjj^fs and tjbh denote the fraction of neutron 
stars and black holes one obtains a contribution {pb/p)NS+BH = {vns + Vbh) ■ 0.3. 
This is still a moderate contribution, but the number of these objects becomes relevant. 
Neutron stars are supposed to provide a portion of less than 1% of all stars, one expects 
a contribution not larger than 10~^. Unless there is an extremely large number of black 
holes the latter will also only give a small contribution. 

Finally let us look at the implications of voids, surrounded by bubble walls. The 
average density of matter is p = 2 ■ 10~^^g cm~^, for the radius of the voids we take 
R = 20Mpc = 6 ■ lO^^cm. From (6.20) we then obtain {pb/p)BW = 2 ■ 10^^. The 
smallness of this contribution is surprising at first sight, therefore it deserves a comment. 
The result can be easily understood. The distribution of matter in the bubble wall is 
homogeneous in the two tangential directions. Only in radial direction it is concentrated 
in a small shell of thickness 2L. The situation is therefore essentially one-dimensional. 
But while in three or two dimensions the Green function of the Laplacian is ~ 1/r and 
~ Inr respectively, i.e. singular at the origin, in one dimension it is ~ r, which is finite. 
Therefore one can perform the limit L — > in the latter case. No problems arise if 
one squeezes the matter in the surface r = i? by taking p{r) = m5{r — i?)/47ri?^. The 
solution hQQ{r) in (6.19) is finite at r = i?, it just has a kink there. The absence of any 
singular behavior for small L for matter concentrated in walls explains why the effect is 
so small. 

We are aware that our result is in striking contrast to the statements of Wiltshire [TS] 
who claims that the presence of large voids should lead to a considerable backreaction. 
We don't feel in a position to comment on this work. But we doubt whether a separate 
treatment of the metric within and outside the voids is appropriate or even legitimate. 

The relative unimportance of backreaction is, of course, by no means a new result. 
It was already found (in a non covariant calculation) by Nelson [H] in 1972, also by 
Wetterich [lOj who's model we essentially used, as well as by many others. But certainly 
the discussion on the (un)importance of back reaction will not end in the near future. 

7 Summary and conclusions 

If one wants to go beyond the simple static approximation which we used for the cal- 
culation of backreaction one is faced with the problem that the metric is not known 
from the beginning. It has to be determined e.g. from a perturbative solution of the 
field equations which couple metric and matter. This calculation can be done in any 
reasonable gauge, e.g. in the harmonic gauge. Subsequently one can transform to our 
optimal gauge, thus removing all unphysical gauge modes. 

The covariant fitting procedure presented here is, of course, not unique because one 
could modify the minimization steps which define the optimization. But it is not at all 
trivial to formulate conditions which do not involve the initially unknown scale factor. 
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can be simply treated perturbatively, and lead to a maximal fixing of the coordinates. 
One could e.g. exchange the order of the steps. We did not perform a systematic 
investigation of all possibilities but chose an order which was motivated by technical 
simplicity. If one would perform steps 2 and 3 at the beginning one would obtain the 
synchronous gauge. The E-L equations for the former step 1 would subsequently become 
rather ugly. There is now no reason for (^gS) and (^G^™) to vanish, one could obtain 
a (small) backreaction already in first order. 

It is worthwhile to emphasize that familiar gauge fixing prescriptions in the literature 
do not fix the gauge in the maximally possible way as it was achieved here. This is well 
known for the synchronous gauge, but it is also true for the Newtonian gauge, even if 
one restricts to time independent transformations. As an example consider the nonlinear 
transformation = 0, — a^r'^ — 2x^ J2k ci^x^-, with a constant vector. This leaves 
Qqq and g^Q unchanged, and only changes A in the decomposition (3.10) to ^4— 4 a^x^. 
Therefore demanding e.g. B = Cm = = is not sufficient for fixing the gauge. 

We came to the conclusion that dark energy cannot be mimicked from backreaction. 
The only quantities which agree are the signs > and < 0. But the relative 
importance of backreaction in the present day universe turned out to be of order 10~^ to 
at most 10~^. Somewhat different parameters than the ones used above, or more realistic 
density distributions inside the sources, would not change the numbers considerably. We 
also gave arguments why a perturbative expansion makes sense, and why neither the 
presence of large voids nor retardation effects should seriously modify the results. Only 
if one widely gives up the cosmological principle one may evade these conclusions, but 
then it becomes hard to derive even semiquantitative statements. 

In order to mimic dark energy from backreaction the small ratio Pb/p should somehow 
increase to the observed ratio Pde/T> ^ 7/3 between dark energy and matter. It appears 
miraculous how this could happen. An even greater miracle would be needed to change 
the ratio Wb = —1/21 to the observed wde ~ —1 for dark energy. 

Acknowledgement: I thank E. Thommes for valuable discussions and for reading 
the manuscript. 
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